Renormalized perturbation theory and scaling for an impurity Anderson model 
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We demonstrate the effectiveness of a generalized renormalized perturbational approach to cal- 
culate the induced magnetization for the single impurity Anderson model with a strong on-site 
interaction, using flow equations for renormalized parameters to scale from a weak correlation to a 
strong correlation regime. We show that, using simple approximation schemes in different parameter 
regimes, remarkably accurate results can be obtained for all magnetic field values by comparing the 
results with those from direct numerical renormalization group and Bethe ansatz calculations. 



I. INTRODUCTION H. 



G a {w,n) 



The impurity Anderson model in the strong interaction 
(Kondo) regime has served as a classic testing ground for 
the development of many-body techniques for application 
to a wider class of models for strongly correlated electron 
systems. This is because exact or very accurate meth- 
ods have led to a comprehensive understanding of this 
mode! 1 ; notably via the numerical renormalization group 
(NRG) 2,3 , Bethe ansatz method 4 ^, and for degenerate 
models the 1/N— and slave boson methods 7,8 . Though 
a direct perturbation expansion in powers of the on-site 
interaction U works well for low and intermediate inter- 
action strengths for the symmetric Anderson modeler— 
but it breaks down in the interesting strong interaction 
regime where the model displays all the strong corre- 
lation physics of the Kondo model (it can be mapped 
into this model via a Schrieffer- Wolff transformation 13 ) . 
The assumption has been, therefore, that the strong cor- 
relation low energy physics can only be obtained using 
a non-perturbative method. However, the low energy 
fixed point corresponds to a Fermi liquid 14 , which can 
be described in terms of quasiparticles, with the single 
particle excitations in 1-1 correspondence of the non- 
interacting system, together with a limited number of 
parameters to describe the inter-quasiparticle interac- 
tions. As a consequence the perturbation theory can 
be reorganized as a renormalized perturbation expan- 
sion, similar to that used in quantum electrodynamics, 
with parameters corresponding to renormalized versions 
of those of the original model. The renormalized pertur- 
bation theory (RPT) provides asymptotically exact re- 
sults for the thermodyamics and dynamics in the low en- 
ergy regim o 15 ' 16 . However, the renormalized parameters 
have to be calculated and in the strong correlation regime 
this would seem to require a non-perturbative technique 
such as the NRG 1 - 7 -. 

Here we demonstrate that this is not the case. We in- 
troduce a generalized scaling scheme, such that we can 
follow a set of flow equations for the renormalized param- 
eters from a regime of weak to strong electron correlation. 
Once these parameters have been determined, there are 
exact relations from which the magnetization, spin and 
charge susceptibilities 1 ^ 1 ^— can be calculated. G a (w,fi) = 

We consider the Anderson model 22 in a magnetic field 



HAM = ^ € d,a d lr d cr + u n d,T n <U W 
cr 

e k,a c k,cr c k,cn 



where ea, CT = £ d — a = ±1 (t,4-), is the energy 

of the localized level at an impurity site, U the inter- 
action at this local site, and T4 j[T the hybridization ma- 
trix element to a band of conduction electrons of spin a 
with energy tk.o — o~g c fiBH/2, where g c is the g-factor 
for the conduction electrons. The resonance width factor 
A ff (w) = 7r^ fe |Vfe !<7 | 2 (5(ai — £k,cr) we take as a constant 
A, corresponding to a flat wide conduction band, and we 
use the notation h = g^H /2. The model can be charac- 
terized by the set of parameters (h, e^, A-j-, Aj,, U), which 
we can regard as a defining a point fi in a configuration 
space. 

The Fourier transform of the retarded single-particle 
impurity Green's function has the form, 



1 



(2) 



where fi) is the self-energy, which we assume to be 

non-singular at lo = so that the low energy fixed point 
of the system corresponds to a Fermi liquid. The renor- 
malized parameters 2 ^, ed, „•(/•*) and A CT (/Lt), which charac- 
terize the quasiparticles of the Fermi liquid are defined 
by 

h\,Ap) = z<r(ed,<7 + E ff (0,/x)), A a (fi) = z a A, (3) 

where z a is given by z a = 1/(1 — 5^(0, ju)). It is useful to 
introduce a renormalized field h(fi) and average effective 
level €d(n) defined by 

g d(M) = 1^2 ^O): ~ h (v) = -\Yl < rl &Al J ')- ( 4 ) 



We can then define a quasiparticle Green's function via 
G (7 (w,/Lt) = z a G a (u, fi) such that 

1 



uj - €d(n) + ah(fi) + iA a (fi) - S CT (w, fi) 



(5) 
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where T, a (uj,fi) is the renormahzed self-energy. The 
effective interaction between the quasiparticles U(fi) 
can be defined in terms of the local total 4-vertex 
Y^){wi,W2, ^3,^4 : fi) at zero frequency, 



1700 = 2^4.^(0,0,0,0 



(G) 



The Lagrangian corresponding to the bare Anderson 
£am(m) can be rewritten in the form, 



£am(a0 = £am(m) + £ ct (A), 



(7) 



where A are the coefficients in the remaining or coun- 
terterm part of the Lagrangian. They can be determined 
from the renormalization conditions, 

E(0,m)=0, E'(0,/i) = 0, f t4 (0,0,0,0:/x) = C/(/i). 

(8) 

These ensure that there is no overcounting of renormal- 
ization effects which have already been included by using 
the renormahzed parameters instead of the bare ones. 

Exact results can be derived for a number of physi- 
cal quantities in the Fermi liquid regime in terms of the 
renormahzed parameters 15 . Calculations can be carried 
out for the spectral functions and behavior beyond the 
Fermi liquid regime using renormahzed perturbation the- 
ory (RPT) in terms of the renormahzed parameters as no 
approximation is made in setting up the expansio n 16 ' 18 . 

The renormahzed parameters can be regarded as defin- 
ing an alternative specification of the model, ft = 
(h(fi),e d (fi), A t (/x), A±(fi),U(fi)). Here we introduce 
the idea of deriving a scaling relation linking the renor- 
mahzed parameters for a model specified by the param- 
eter set ii\ to one specified by the parameter set fi2, 
/Xi(/Lti) — > /X2(a*2)- We consider a renormahzed pertur- 
bation expansion for G a (uj, fi 2 ), in which we use the non- 
interacting quasiparticle propagator for the system with 
fi = fix instead of that for the system with (1 = fi 2 - 
Expressing this Green's function in the form, 



G a (u},fi 2 ) 



2<r(Ml) 



(g£ 0) (u;,mi)) -1 - S CT (w,^i,^2 - Ml) 



(9) 

effectively defines the self-energy S CT (u;, fix, fi 2 — fii). The 
corresponding Larangian takes the form, 

£am(/-*2) = £am(Ai) + £ ct (Ai) + C cx (fi 2 - Mi), (10) 

where there is an additional term £ ex (/2 2 — Mi) due to the 
difference between the models with different parameter 
sets. The counterterms in £ ct (Ai) are required to satisfy 
the conditions (|5]) for ft = fix. 

By equating the inverse of G CT (o;,M2) from Eqn. <j9j) 
with that derived from Eqn. ([5]) and its derivative at 
uj = 0, we can relate the renormahzed parameters at 
Hi to those for fi2- For fix = fi, fi 2 = M + <^M> these 
equations, take the form, 



where z a (fi,Sft) = (1 — £^(0, fi, Sfi)) 1 , and 

e d (fi+5fi)-ah(fi+Sfi) = z a (e d (ft)-ah(fi)+'S a (0, fi, 6fi)). 

(12) 

The self-energy S CT (a;, fi,Sfi) and its derivative are zero 
at u = for Sfi = , so we take them to low- 
est order in Sfi, which in general will be first order in 
8fi. We write them as S CT (0, fi, Sfi) — ax 1<T (fi,Sfi) and 

S^.(0,Mj^m) — Q! 2,cr(Mi ^M)- Working to lowest order in 
5fi we can separate (|12p into two equations, 

e d (fi + Sfi)-e d (fi) = (13) 
e d (fi)a 2 , + (fi,Sfi) - h(fi)a 2 -(fi,Sfi) + ax,+(fi,Sfi). 



h(fi + Sfi) - h(fi) 
h(fi)a 2 ,+ (fi,Sfi) 

where 



e d (fi)a2-(fi,Sfi) 



(14) 

ax,-(fi,Sfi). 



OH,±{l*, <*M) = 0-5(ot,t(M) ^M) ± "x,4-(M) <^M))- ( 15 ) 

To apply these scaling equations to calculate the renor- 
mahzed parameters in a strong correlation regime two 
conditions need to be satisfied: (i) we need a known set of 
renormahzed parameters corresponding to fix, Mi ~~ * Mii 
and (ii) an accurate method of calculating the renormal- 
ized self-energy E CT (w, fi, 5fi) in the low frequency range, 
ie. its value and derivative as u> — > 0. The results should 
be independent of the particular trajectory chosen for 
generating fi2 from fix- Note that in contrast to the 
Wilson renormalization approach, where states are in- 
tegrated out, one can define an inverse transformation 
M2(M2) Mi(Mi)- Condition (i) can be satisfied by 
choosing fix in a weakly correlated regime where per- 
turbation theory can be applied. The more difficult step 
is to find an approximation to satisfy condition (ii). 

In particular cases one can use additional information 
derived from relevant Ward identities. We consider the 
case where we change only the parameter fix — h. Using 
the Friedel sum rule 2 ^ and a Ward identit y 10 ' 12 we derive 
the result, 

a lt<T (h,Sh) = -<T(l + U{h)p- a (0,h))5h+O((5h) 2 ), (16) 



where 



A(m) 



*{u:-e d (fi)~h{fiy + &(fiy 



(17) 



The details of the derivation are given in the Appendix 
A. We can use this result in the flow equations. For 
h = 0, e d (Sh) and a 2:a (0,0, Sh) are of order (6h) 2 , and 
from the scaling equations we find h(h)/h — > R where 
R = 1 + U(0)p(0, 0), is the Wilson ratio, which is an 
exact result. 

We can derive similar relations in the case where we 
change fi 2 = e d , 



A a (fi + Sfi) = z a (ji, Sfi)A a (fi), 



(11) 



ai^{e d ,8e d ) = (1 - U(e d )p^ a (0, e d ))Se d + 0((Se d ) 



(18) 
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where again details of the derivation are given in the 
Appendix A. If we change both and h, then these two 
results can be combined in the scaling equation to first 
order in both Sed and Sh. We use these results in the 
scaling equations to show that suitable approximations 
can be derived to satisfy both conditions (i) and (ii) in 
particular parameter regimes. 



II. SCALING WITH MAGNETIC FIELD 

We show results first of all for the particle-hole sym- 
metric model. In this case we can satisfy condition (i) 
by considering the system in a very large magnetic field 
so as to suppress the spin fluctuations that lead to the 
strong correlation effects in this regime. We then test the 
accuracy of the approximation used in step (ii) by calcu- 
lating the induced m(h) magnetization at T = 0, which 
is given in terms of the renormalized parameters £d,<r and 
A<r by 



m(h) — a tan 1 ( -i- 

K ' 2tt 4- V A„ 



(19) 



and compare the results with the corresponding NRG 
and Bethe ansatz results. For particle-hole symmetry 
the formula for the magnetization given in Eqn. (|19[) 
simplifies as Id.a = —oh and A a is independent of o. 



FIG. 1: (Color online) A diagram for the (renormalized) self- 
energy involving a repeated scattering in the particle-hole 
channel, corresponding to scattering with a spin flip. In the 
initial calculation in very large magnetic fields, the dashed 
line represents the 'bare' interaction U and the full lines are 
the propagators calculated in mean field theory in an applied 
magnetic field h. In the renormalized perturbation theory the 
interaction parameter is Up h and the propagators are for the 
quasiparticles in an effective field h and include the mean field 
insertions. 



Accurate initial renormalized parameters in the ex- 
treme large field regime h > U can be calculated simply 
using the original Anderson mean field theory 22 . These 
can then be improved upon using the set of RPA spin 
flip scattering terms to calculate the self-energy and the 
corresponding value of U. This is illustrated in Fig. Q] 
and Fig. [2] where the propagators correspond to those 
for the "bare" electrons with mean field insertions and 
the interaction parameter represented by a dashed line is 
the bare interaction U. 

To extend the results to lower magnetic field values we 
use the scaling approach described in the previous section 















u = 


! + 




+ : 


i 
i 


+ 








1 

^3>— 


i 





FIG. 2: (Color online) The diagrams, at zero external fre- 
quency uj = 0, used to calculate the effective interaction ver- 
tex U in the particle-hole symmetric regime in an applied 
magnetic field. In the inital calculation the propagators cor- 
respond to mean field theory and the interaction parameter 
(dashed line) to U. In the RPT calculations, the propagators 
correspond to free quasiparticles and the interaction parame- 
ter (dashed line) is U^ h . 



and reduce the applied field to h — Sh but use the renor- 
malized parameter at the field value h. We then need 
a suitable approximation to calculate the renormalized 
self-energy tl a {uj,h,Sh). We have shown that accurate 
results can be obtained for the renormalized parameters 
by using the same set of diagrams for the self-energy 
and the vertex U as in Fig. [T] and Fig. fJl where now 
the propagators correspond to the free quasiparticles 24 . 
The counterterms have to be taken into account and as 
a consequence the interaction parameter is modified to 
U — A3, which we denote by f/* ft , due to the countert- 

erm A3. Given the starting values of h(h) and A(/i) (for 
particle-hole symmetry = 0) the corresponding value 
of Up h (h) can be deduced from 



m(h) 



= —tan 

7TAI 



h{h) 



n ph (o,/i) 



ui h % h (Q,hy 



(20) 



where YLp^iw, h) is the transverse dynamic spin suscep- 
tibility due to the non-interacting quasiparticles 16 . The 
full expression for Hphiu), h) is given in the Appendix B, 
and its value at u — is given by 



n, 



lph(0,h) = -Ltan 
7m 



irh 



\A(h)J 



(21) 



Eqns. (|2TJ)) and (|2"Tj) imply that h(h) takes a generalized 
mean field form h(h) — h + U^ h (h)m{h). 

The new parameters h and A are then calculated at 
the reduced field h — Sh from the scaling equations (fTTj) 
and (fPo)) using Eqn. (fT5|) with U(h) given by 



U(h) 



TP 



1 



m ph iL ph (o,hy 



(22) 



corresponding to the diagrams in Fig. [2l Once h and A 
at h — Sh have been calculated the corresponding values 
of Up h and U can be deduced from Eqns. (f2TI)) and (|2"2")l 
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and the process repeated to reduce the magnetic field in 
a sequence of steps to zero. 

Results for the induced magnetization for the case, 
U/ttA — 3, 7rA = 0.1 and a conduction band halfwidth 
D = 1, are shown in Figs. [3]and[H together with those 
from a direct NRG calculation. The calculated Wilson 
ratio, R = l + U(0)/irA(0) = 1.99, showing that to an ex- 
cellent approximation U(0) = 7rA(0) so that there is only 
a single energy scale corresponding to the Kondo tem- 
perature X"k- The Kondo temperature Tk is defined in 
terms of the zero field spin susceptibility \s a t T — via 
Xs = ( gm ) 2 S(S+l)/3T K with S = 1/2. The value of T K 
is 0.00800 which is very close to the Bethe ansatz result 
for these parameters 0.00 8 05 25 . The results in Fig. [3] are 
plotted as a function of ln(/i/Tk) over the whole range of 
magnetic field values to the saturated value m(h) 1/2. 
It can be seen that there is excellent agreement between 
the results deduced from the RPT scaling equations and 
the ones determined from a direct NRG calculation^. In 
Fig. Q] we plot and compare the NRG results as a func- 
tion of /i/Tk up to a value h = 4Tk) and also with Bethe 
ansatz results^ for the s-d model over this range. It 
can be seen again that on this scale there is still excel- 
lent agreement between the results. For larger values of 
the magnetic field h > 8Xk for U/irA = 3 the approach 
to saturation is different for the s-d and Anderson mod- 
els because in this range the impurity charge excitations 
play a role for the Anderson model which are not present 
in the s-d model^. 

In earlier results 2 ^ using this approach we have given 
more extensive comparisons with Bethe ansatz and NRG 
results including calculations of the self-energy and dy- 
namical spin susceptibilities. We note a scaling approach 
using a magnetic field using the functional renormaliza- 
tion group method (fRG) has recently been developed 
by Streib, Isidori and Kopietz 26 for the strong correla- 
tion regime for the particle-hole symmetric model. 




FIG. 3: (Color online) The magnetization m(h) as a func- 
tion of a function of fn(/t/Tk) calculated from the RG flow 
in the RPT (dashed line, circles) for the symmetric model 
with U/ttA = 3, TV A = 0, 1 and T K = 0.008 compared with 
the corresponding results of a direct NRG calculation (dashed 
line). 




FIG. 4: (Color online) The magnetization m(h) as a function 
of h/T K calculated from the RG flow in the RPT (full line, 
circles) for the same parameter set as in Fig. [3]compared with 
the corresponding results of a direct NRG calculation (dashed 
line) and Bethe ansatz results^ 1 - 5 (dot-dashed, crosses) for the 
s-d model. 



III. SCALING RELATIVE TO THE CHEMICAL 
POTENTIAL 



Well away from particle-hole symmetry with either a 
low or high density of particles or holes, we have another 
regime where the correlation effects are weak. Here we 
can also apply mean field theory to calculate the initial 
values for the renormalized parameters, and then extend 
the calculation of the self-energy using standard prtur- 
bation theory in U to take into account the diagrams 
with repeated particle-particle scattering (see Fig. [SJ. 
We consider this situation initally without an applied 
magnetic field and introduce a magnetic field later. This 
regime was first studied by Schrieffer and Mattis2£. They 
calculated an effective interaction f7 G ff from an expansion 
in powers of U corresponding to the diagrams shown in 
Fig. 1 



U cS = 



u 



1 + -^tan" 



(23) 



where is the renormalized level calculated in mean 
field theory. The value of U e g can be identified as U and 
is in good agreement with the values deduced from the 
NRG in the regime well away from particle-hole symme- 
try. To scale from this regime to a more strongly corre- 
lated one we decrease the value of | j , which is equiva- 
lent to changing td relative to the chemical potential. We 
use, as an approximation for the renormalized self-energy 
in the scaling equations for part (ii) of the calculation, 
the same set of diagrams with repeated particle-particle 
scattering in the RPT but using quasiparticle propaga- 
tors with a renormalized interaction in this channel U m . 
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The initial value of U — A3 

1 - n d ^ - n d ^ _ 
2E d 



U pp is calculated from 



1 



tan 

nE d \A 
tan" 1 (I) 



7re d 



-(*)' 



(24) 



and the initial value of U from (|23j) . We have introduced 
the notation E d = e d + U/2, so that E d = corresponds 
to the model with particle-hole symmetry. 

Once the values of e d and A have been calculated from 
the scaling equations at e d — Se d , new values of U pp can 



be calculated from Eqn. (|24p 
values of U deduced from 



at e d — 5e d and the new 



FIG. 5: (Color online) A diagram for the (renormalized) self- 
energy involving a repeated scattering in the particle-particle 
channel. In the initial calculation in the low density regime 
of either particles or holes the dashed line representing an 
interaction vertex is the 'bare' interaction U and the full lines 
are the propagators calculated in mean field theory. In the 
renormalized perturbation theory the interaction parameter 
is Upp and the propagators are for the quasiparticles in an 
effective field id and includes the mean field insertions. 



u 



Upp 



^tan 



■'(*)' 



(25) 



corresponding to the diagrams in Fig. [BJ Note in this 
case, in scaling with e d in the flow equations, we are 
changing the model itself rather than an external field. 

To test the range of validity of this approximation 
we calculate the renormalized parameters from the flow 
equations and compare them the corresponding values 
deduced from an NRG analysis of the low energy fixed 
point for the case U/ttA = 3, 7rA = 0.1 over a range 
of values of e d . The approximation based on this set of 
diagrams works well in the range \E d \ > U/2. The re- 
sults are shown in Fig. [7] Corresponding results for the 
static spin susceptibility due to the impurity, which can 
be evaluated from the renormalized parameters using the 
exact result, 



u = 











! + 




+ ! 


1 
1 








1 



+ 



FIG. 6: (Color online) The diagrams, at zero external fre- 
quency uj = 0, used to calculate the effective interaction ver- 
tex U in the regime \E&\ > U/2. In the initial calculation the 
propagators correspond to mean field theory and the inter- 
action parameter (dashed line) is U. In the RPT calculation 
the propagators correspond to free quasiparticles, include the 
mean field insertions, and the interaction parameter (dashed 
line) is U pp . 



Xs = 



(5MB ) 2 - 



p(0,e d )(l + U(e d )p(0,e d )), (26) 



are shown in Fig. [5] The values of \s are shown on 
a logarithmic scale and are in good agreement with the 
NRG results. 

To calculate the induced magnetization m[h) in a mag- 
netic field in this regime the formulae become more com- 
plicated. The equation for U pp generalizes to 



n d ,f + n d ,i 
2E d 



1 



n pp (o) 



Upp^pp (o) 



and for U, 
where 

n pp (o) 



u = 



u, 



I'll 



I - UppTlpp(0) 



(27) 



(28) 



2i d 



7T4e- 2 d + (A t -A ; )2^ 



^ tan -1 



A„ 



27T4^+(A t -A 4 ) 2 



In 



5i,t 



_ e d4 



(29) 



The flow equations are solved as a function of e d for 
a given magnetic field value h, and the renormalized pa- 
rameters used to deduce the induced magnetization using 
Eqn. (|T5|) . The results for m(h) as a function of \n.[h/T^) 
for the case U/ttA = 3, ttA = 0.1 and E d = -0.15 
are shown in Fig |9] This is a significantly correlated 
regime with U/U = 0.29 and Up(0) = 0.54. the value of 
Tk = 0.052 and is defined in terms of the T = suscep- 
tibility as earlier. A comparison is made with the corre- 
sponding results from a direct NRG calculation. The two 
sets of results can be seen to be in excellent agreement 
over the full range of magnetic field values. In Fig. [TUlwe 
compare the two sets of results for smaller values of the 
magnetic field, over a more physically accessible range, 
as a function of h/T^. 

As already noted, the approximation used for the cal- 
culation of the renormalized parameters, based on the 
particle-particle scattering diagrams, breaks down at the 
point I .Ed I ~ U/2. The difficulty in extending the calcu- 
lation beyond this point to the particle-hole symmetric 
point I Ed I = is that, in the strong correlation regime 
U/ttA 3> f, the effective interaction in the particle- 
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FIG. 7: (Color online) The renormalized parameters id/nA, 
A/ A and U /-kA as a function of Ed calculated from the RG 
flow in the RPT for the model with U/tvA = 3, ttA = 0.1, 
compared with the corresponding results of an NRG calcula- 
tion. 



should go to zero as the magnetic field is reduced, How- 
ever, in the numerical calculations a small difference per- 
sists as h — » 0, so an improved or alternative approxima- 
tion is required to obtain completely satisfactory results 
in this parameter regime. This is currently being inves- 
tigated. 
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FIG. 8: (Color online) The static susceptibility (on a log scale) 
in the absence of a magnetic field as a function of Ed cal- 
culated from the RG flow in the RPT for the model with 
U/ttA — 3, n A — 0.1 compared with the corresponding re- 
sults of a direct NRG calculation. 



FIG. 9: (Color online) The magnetization m(h) as a function 
of ln(h/T K ) calculated from the RG flow in the RPT (full line, 
circles) for the model with U/ttA = 3, ttA = 0, 1, E d = -0.15 
and Tk = 0.052, compared with the corresponding results of 
a direct NRG calculation (dashed line) 
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particle channel U pp becomes very large so as to sup- 
press the charge fluctuations. In the strong correlation 
regime near half-filling U p(0) — > 1 so that the impurity 
charge susceptibility, which is proportional to 1 — Up(Q), 
is effectively zero. It can be seen from the scaling Eqn. 
(O and Eqn. (JTHJ) that when Up(0) = 1 the first or- 
der term in 6ed is zero. Near particle-hole symmetry for 
strong correlation, the main effect is simply a change in 
the Kondo temperature Tk, as Up(0) — 1 is maintained 
in this regime, and Ty^(Ed) — Tk(0) oc E^. 

An alternative strategy to calculate the renormal- 
ized parameters in the strong coupling regime close to 
particle-hole symmetry would be to scale from large to 
small magnetic fields, based on the particle-hole scatter- 
ing diagrams in Fig. [T] and Fig. [2] Away from particle- 
hole symmetry with moderate to large magnetic fields 
there is a difference between Af(h) and Aj,(/i) which 



FIG. 10: (Color online) The magnetization m(h) as a function 
of h/T K calculated from the RG flow in the RPT (full line, 
circles) for the same parameter set as in Fig. [9]compared with 
the corresponding results of a direct NRG calculation (dashed 
line) 



IV. CONCLUSIONS 

We have shown how, based on a relatively simple set 
of diagrams within a renormalized perturbation theory 
(RPT), accurate results can be obtained for the magne- 
tization and susceptibilities in different parameter regions 
of the impurity Anderson model with strong electron cor- 
relation. The RPT approach is a physically transparent 
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one where we follow the flow equations for the renormal- 
ization of the quasiparticles from their almost bare val- 
ues in a weak correlation regime to strongly renormalized 
ones in the strong correlation regime. In contrast to other 
renormalization techniques which involve eliminating or 
integrating out of higher energy states, such as the NRG, 
we need to consider the flow of a restricted number of pa- 
rameters only, those required to specify the low energy 
behavior, such as A and h, or A and e d . This approach 
opens up the possibility of applications to a much wider 
class of models involving strong electron correlation. 
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V. APPENDIX A 

We consider the derivation of the result in Eqn. (fTBj) 
where we change only the parameter p\ = h. We can use 
the Friedel sum rule, 

nd,a{h + Sh) = - 

_I tan -i ^ e d -(T(h + 5h) + Z a (0,h + 6h) ^ (3Q) 

Expanding this to first order in Sh, 
nd,a(h + Sh) = n d , a {h) 

+a8h Pa (0, h) (l - o- d^f^ + 0((Sh) 2 ), (31) 



Oh 



where 



1 



ir (w - e d + ah - S CT (0, h)) 2 + A 2 ' 
We also have a Ward identity in the form, 



(32) 



=o-<? ETT = P-cr{0,h)r t ^(0, 0,0,0), 



ouj ah 

which can be rewritten in the form 
dZ a (0,h) 



(33) 



1 - (7- 



^(0,^^(0, 0,0,0) + —. (34) 
oh z a 



Using this result in Eqn. (|4Tfl) we obtain 

nd, a (h + Sh) = nd,a{h) 

+a8hp a {0, h)(l + U(h)p- a (p, h)) + 0({Sh) 2 ),{35) 

where we have used the fact the p a (0, h) = p CT (0, h)/z a . 
We also have 

n d ,a(h + 8h) = i 

_ 1 x ( i d {h) ~ <J~h(h) + t a {0,h,8h) \ 
^ I A a (h) ' 



Expanding this to first order in Sh, 

nd : a{h + Sh) = n d ,a{h) 

-aShp a (0, h)£ a (0, h, Sh) + 0((Sh) 2 ), (37) 

Equating the two expressions we find 

a\, a (h,8K) — Yi a {0,h,5h) = 

-cr(l + U(h)p-*(0, h))5h + 0((Sh) 2 ). (38) 

Using a very similar argument we can derive the result 
in Eqn. (|18j) where we consider a change 5p2 = 



/ lt n 1 1 -l {ed + 8ed + Za(0,ed + 8e d ) 

n d ,a{ed+d€d) = tan 

Z 7r \ Za 



(39) 



Expanding this to first order in Sed, 
n d ^{e d + Se d ) = n d ^{e d ) 

-8e dP a(0, e d ) (l + ^g^ ) + 0((Se d ) 2 ). (40) 

We also have a Ward identity in the form, 
ggoKfr) , , dE a (0,e d 

n \U — 0"T 

ouj oe d 

which can be rewritten in the form, 
d^(0,e d ) 



/ 9_ <7 (o,/i)r ff ,_ CT (o, o,o,o), 

(41) 



1 



1 + " " v "'""" = -P-A0, e d )r CT ,_ ff (0, 0, 0, 0) + — . (42) 
Substituting into Eqn. (|4*U|) we find 

nd,<jy{td + Sed) = rid,a(ed) 

-5e d p*{0, e d )(l - U(e d )p- a {0, e d )) + 0((Se d ) 2 ){^) 

We also have 

/ . r n 1 1, -x / ^d + ^a(0,ed,8ed) \ 

n d A e d + Se d ) = tan =-— . 

2 7T V A(e d ) / 



Expanding this to first order in Sed, 



(44) 



n dirT (ed+Se d ) = nd. < y{e d )-p CJ {0, e d )t, a (0, e d , 8e d )+0((Se d ) 2 ), 

(45) 

Equating the two expressions we find 

ai^(e d , Se d ) = £<t(0, e d , 8ed) = 

(1 - U(e d )p-*(0, e d ))Se d + 0((5e d ) 2 ). (46) 

VI. APPENDIX B 

For the model with particle-hole symmetry the dy- 
namic transverse spin susceptibility for the free quasi- 
particles (particle with spin j" and the hole with spin \) 
in a magnetic field is given by 



r 7T 



(47) 
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for cj > 0, where 
F ph (cu) 

K ph (oj) = 



1 



^ln 



cu + 2h 



oj + h + iA 



1 



-^ln 



uj + 2h + 2iA 



-h + iA 



uj + h + iA 



h + iA 



(48) 



(49) 



The causal propagator U p h(uj) is an even function of uj 
in the absence of a magnetic field. 

The two-particle propagator for free quasiparticles 
fl pp (uj) is given by 

n PP H = ±-(f™„ + fz,„ + K% + kz,„) (50) 



for uj > 0, where 



1 



uj - 2l d + iA a - iA- a 



-In 



and 



la.-e + iA- a 
(51) 



1 



uj - 2e d + iA a + iA- a 



-In 



iA n 



-l d ,o + iA 



(52) 

In the absence of a magnetic field the causal propagator 
IIpp(w) is an even function of uj. 
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